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ABSTRACT 



Numerical weather analysis relies heavily on numerically computed 
winds both for direct output and internally for calculating advection and 
other purposes. It is the intent of this paper to compute, analyze, and 
compare the geostrophic, the non-divergent, the divergent, and a com- 
bination of the latter two. Where there occurs a significant difference 
between these numerically derived winds, an attempt is made to com- 
pare each to the observed winds. This study is made for all standard 
levels through 300 mb. Data and several computer programs were fur- 
nished by the U. S. Navy Fleet Numerical Weather Facility, Monterey, 
California . 

The writer wishes to express his appreciation to Professor George 
J. Haltiner of the U. S. Naval Postgraduate School for his guidance, 
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1 . Introduction 



One of the primary products of any numerical weather facility is 
wind data for operational use and for further numerical computations 
such as advection. It is the purpose of this paper to analyze and 
compare several types of numerically derived winds. The types con- 
sidered are the geostrophic, Vg, the non-divergent, Ve, the divergent, 
Vx, and finally the total wind, Vt, which is a combination of the di- 
vergent and non-divergent parts of the wind. These winds are cal- 
culated by solution of appropriate equations using data from the U.S. 
Naval Fleet Numerical Weather Facility, Monterey, California (FNWF). 
Computations were made on the U. S. Naval Postgraduate School CDC 
1604 digital computer. 

Comparison of these winds is made at all standard levels through 
the 300-mb level by calculation of speed and geographical direction on 
separate maps for each level covering the Northern Hemisphere. Dif- 
ference fields of direction and speed are made between each wind and 
the geostrophic wind. Vertical cross sections were computed and 
studied for each of the winds and for the difference fields. Further com- 
parisons and analyses are made by using these three winds to advect 
absolute vorticity. Advection of temperature was also computed using 
Vx. 
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2. Background 



The geostrophic wind is computed from the well known formula 




( 1 ) 



utilizing a standard program of FNWF, This program generates fields 
of geostrophic wind directions in degrees and speed in knots. All com- 
putations were done on the FNWF square grid centered on the North 
Pole and consisting of 3968 with a grid interval of 381 km at SON. 

In this program the sine of latitude, as required in computation of the 
Coriolis parameter, f, was set equal sine 15 degrees. The latter points 
were not utilized in the evaluation of results. 

Geostrophic wind , Vg, is utilized as a basis of comparison for 
later derived non-divergent, Ve, and total, Vt, winds. 

The Helmholtz theorem allows expression of the wind as the sum 
of the divergent and non=-divergent components. 



Winds of this form will be ••‘ttiii referred to as total wind, Vt.. 

The potential function X for the divergent part of the wind is obtain- 
ed through relaxation of the Poisson-type continuity equation 



XV>f +■ VX. 



( 2 ) 




(3) 
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The stream function, 't, for the non-divergent portion of the wind is 
approximated by solution of the balance equation 



( 4 ) 



Here, <S> is the geopotential. 

The vertical component of vorticity, C > is the Laplacian of the 
stream function 



f = 7^ r 



and absolute vorticity, 77 , is the sum of fand the vertical component 
of the earths vorticity, f, 

7j= f . ( 6 ) 

Vertical velocity, co, is computed by the diagnostic equation 



y^fcroo)f^ -;^[y^J(OjT)~J(D,V^T) ^ J (T,7}jJ 

developed by Haltiner, Clarke, and Lawniczak [l]. 

For convenience the stream function has been redefined so as to be 
in units of height, 

f ^ ^ E. (8) 

The geostrophic advection of absolute geostrophic vorticity is com- 
puted from the expression 

-f Jcd.'HsJ (9) 
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and the advection of tIq by 

_ 5 






(10) 



The advection of tJq by Vx is computed from the formula 

d 9?e 

Ay ‘ 



\Vi • V ^ 



( 11 ) 



(10) and (11) are combined to represent advection of rjQ by Vt and 
compared to similar advection using Vg and Ve, 
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3 „ Procedures 



The sequence of calculations programmed started with solution of 
the geostrophic wind equation in finite difference form 



Fields of Ug and Vg relative to map projection were generated and convert- 
ed into fields of wind direction in geographical headings at each grid 
point. Likewise, fields of isotachs were produced and printed in con- 
toured map form. These calculations and maps were made for all stand- 
ard levels to 300 mb for this and all subsequent programs. 

The balance equation, (4), as progra mmed by M. Reese of FNWF, 
was then solved and the fields of E values were stored on magnetic 
tape. The iterative scheme of this program does not always converge at 
300 mb or has an excessive running time, both of which limited the data 
produced. The E values were then used as input to the geostrophic 
wind program with the appropriate substitution of f for f. These "balanced" 
winds were computed and stored in the same manner as Vg. 

Divergence fields are taken as input for solution of (3), Here, as 
with X is redefined as 




(12) 




as) 
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so that the velocity potential has units of height and therefore is com- 



patable with E and D foruss-a ii (12) with minor changes. Expressing 



the first term of equation (3) in finite difference form gives 



2 . 

V X ^ 



3 



^ uj 

dp 



(14) 



which may be solved by use of an FNWF Poisson relaxation program. 
The derivitive “S oo/ S P which is the horizontal divergence of the wind 
is evaluated by a finite difference approximation as programmed by 
Alden and Rosenberger L2^ 

Th*v^ fields are then used to compute the divergent wind component 
by the finite difference equation 









(15) 



'X - 

Again the u and v components are stored and maps of direction and of 
isotachs are produced. 

Now the components of Ve and Vx are added and a field of these com- 



bined components produced as follows 
u Ux = ut 



(16. a) 



V y + Vx = vt (1 6 . b) 

As before these component fields are stored and maps of each level print- 



ed. 



Programs to aid comparison and analysis are now written. Analysis 
of directional differences is made by using the direction fields of Vg as 
a basis and subtracting directions of Ve and Vt from it. This proved to 



be a challenging problem for several reasons. First, directions are 
by nature discontinuous at singular points on isogons maps and also 
in ^neral at 360 degrees. Second, direction tends to vary widely 
in the vicinity of small scale features that are not concentric between 
the various fields compared. 

The discontinuity problem was solved successfully with a rather 
intricate program. However, the singular points and non-concentric 
small scale features still caused a large gradient of direction and 
difficult contouring to plot or read. This problem was alleviated by re- 
placing a value of minus zero at all points on the directional difference 
maps at which the wind speed was less than five knots. This improved 
appearance of the maps and facilitated evaluation. Moreover the min- 
us zero values still indicated those areas where singular points had 
occurred. The 5=knot criterion proved to be of no detriment in subsequent 
use. 

Speed-difference fields, IGIE and IGIT, are produced by the same pro- 
gram as directional difference. The computation of speed differences in- 
cluded a statistical analysis and grid maps of output included the pillow, 
P, and root- mean- square difference, RMSD. These terms are defined in 
APPENDIX I. 

Another series of programs were written to calculate the advection of 
absolute vorticity, V , by the various winds previously computed and to 
compare the results. 
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MM. 



First, fields of D values are used as'input to evaluate the geostro- 
phic vorticity -l/^ ^ ^ m 2. 






(17) 



with an FNWF program. Then the resulting fields from (17) are used to 



obta in 






(18) 



\Miich is plotted on standard grid maps. This program is then internal'= 
ly modified to solve 










(19) 



and 






( 20 ) 



Next previously stored components of Vx are used to compute the 
advection by the divergent wind as follows; 




This result is added to (20) to produce fields of TADV» 

Next these three advection schemes are compared using DADV as a 
basis and forming difference fields similar to those used for direct com- 
parison of the winds . 

Typical computer operating time including efficient handling of tapes 
and printing of output on an IBM 1402 printer is shown in Table I. 

TS:'es shown are for one reporting period of one day at all levels. 
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Table I 






Program 


Run time in minutes 


Solution of Balance Equation, (4) 


46 


Wind speed and directions from various inputs 


18 ^ach) 


Solution of equation (14) for y 


30 


Solution of (16. a) and (16. b) 


12 


Comparison of winds 


35 


Advection of absolute vorticity with all winds 


35 


Comparison of advections 


20 


Preparation of grid field for FNWF plotter 


4 (each) 


Total 


200 minutes . 



It is apparent that operating time is a significant factor in restricting the 
amount of data collected. 
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4. Results 



Results and conclusions are based on analyses of data selected 
from the period February 10 to 20, 1964. Only the United States is 
used in figures as it represents an area of reliable data. Using a 
small area for figures also allows display of finer detail. For illustra- 
tion 0000 Z 16 Feb 1964 was selected as it contains a wide range of 
synoptically interesting features over the United States, 

Figure 1 depicts the pressure analysis (dashed lines) superimposed 
on the stream function analysis at 500 mb. The relative height differ- 
ence at the semi-enclosed low center south of the Great Lakes is typ- 
ical of the 500-mb level throughout the period studied. Figure 2 shows 
stream function analysis of the same period. Comparison with pressure 
height analysis i ndicates that the patterns are similarly related at other 
standard levels. 

Wind speeds computed from 500-mb height and stream function analy- 
sis are shown in figure' 4, Dashed isotachs enclose centers of maximum 
and minimum speeds. Speeds for other standard levels are shown in 
figures 2 and 3 where dashed lines are isotachs. Large numerals indicate 
values of minimum and maximum speeds. Analysis of speeds of these 
two winds shows there is a significant difference between them. Dashed 
lines in figure 5 (lower) are isolines of difference in wind speeds, geo- 
strophic minus non-divergent. Similar analysis of other levels are shown 
in figures 6 and 7 . 
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Figure 5 (upper) shows differences in direction for 500 mb, Dg - De, 
Solid lines show the 500-mb pressure height analysis, alternate dashed- 
dot lines are contours of negative 5-kt speeds, and dashed lines are 5=kt 
contours of positive difference. The area marked > 30 is almost com- 
pletely within an area of wind speeds of less than 5 kt. Difference in 
directions are not shown for other levels but are easily located by refer- 
ence to figures 2 and 3, Pillow and RMSD values for the hemispheric 
maps is given in figures 5 through 7. Note that P indicates that on the 
hemispheric scale Ve speed is lower at upper levels, e.g. by a factor 
of approximately 10 per cent at 300 mb. Except in regions of very sharp 
curvature or near singular points, DGDE, figure 5 (upper), is generally 
less than ten degrees at this and all other levels. A more detailed 
analysis was made of gridpoints and actual winds from 185 soundings 
in areas where DGDE was greater than 10 and both Ve and Vg exceed 
5 kt. The graph in figure 8 shows the relationship found. It must be 
pointed out that interpolation within grid points was required to use 
actual soundings and that actual soundings are transmitted to the nearest 
10 degrees only. In an effort to at least partially reduce effects of fric- 
tion and local terrain effects the surfeice winds were not included in the 
data used for construction of this graph. "Va" is the actual wind as taken 
from the soundings and/or streamline analyses. A similar detailed analysis 
of speed differences was made using 210 soundings from areas where 
IGIE exceeded 10 knots. Results are shown in figure 9. 
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Maps of X for 1000 mb (lower) and 300 mb (upper) are shown in figure 
lOo Dashed lines are 5 meter contours, solid lines are 2-kt intervals 
of wind speed, Ix. Figure 11 shows '<and Ix for 500 mb and 850 mb. 
These maps were plotted using an FNWF program designed for another 
purpose so that the contour labels have no intrinsic meaning. Contour 
interval is 10 meters. Note that the pre-trough region in low levels 
is associated with a x maximum (point of Vx inflow) at upper levels and 
a X minimum (point of Vx outflow) as would be reflected in divergence 
and convergence patterns associated with such height patterns, 

Ix can be seen to be maximum at the 300-mb and 1000-mb levels. 
Note that these maximums are of the same magnitude, 8 to 10 knots, 
but that directions are opposite. As might be expected, speeds are min- 
imum at 500 mb as a reflection of the minimum values of divergence 
known to exist there. 

At 300 mb even maximum values represent a very small percent of 
the wind speed, seldom over 10 percent and generally less than three 
percent of Ve. With the lower wind speeds found at lower levels these 
same magnitudes of Ix are found in some instances to equal 45 percent 
of Ve. Such a case is depicted along the Gulf Coast in figure 8 (lower). 
It is of further interest that a deepening, fast-moving cold front was 
located in the same area almost directly under the 1000-mb Ix maximum 
center. This front moved rapidly up the East Coast and into the North 
Atlantic. A center of maximum Ix was associated with this front through- 
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out the traversal of the United States but speed of movement and the 
short period of accurate data prevented detailed study of this associa- 
tion o 

Combining Ve and Vxj Vt was formed. Figure 12 (upper) shows the 
effects of this combination at 850 mb; 12 (lower) depicts 1000 mb. 

Solid lines are isotachs of Vt, dashed of Ve. The effect of this 
combination at middle and upper levels is so small as to be undetectable 
in this type analysis and is not shown. The effects at 1000 mb are quite 
evident. Detailed analysis of these proved that here there is a change 
in location of relative differences but no change in magnitude of the 
difference. Terrain and friction effects preclude any meaningful com- 
parison of the small deviations from Ve to actual wind data. Further 
attempts to contrast Vt with Vg only further emphasizes that there is no 
significant difference between Vt and Ve. 

Figures 12 and 13 depict the results of using Vg and Ve for advection 

of absolute vorticity. Solid lines are lines of zero advection using Vg, 

DADV; dashed lines are for Ve, EADV. The 300-mb level is omitted but 

showed no marked difference between the two advecting winds as it too 

showed practically congruent patterns. The numerals in the ratios shown 

are located at points of local minimum and maximum value, the upper num- 

=>7 ~1 

eral is for DADV. Units are 10 sec . These values of change of 
vorticity are typical of the period investigated and showed DADV to be 
larger at both extremes by a factor of 5 to 10 units. Computing the pillow 
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and RMSD for all the standard levels gave a pillow O in every case and 
RMSD of 2 to 7 units with largest being at 300 mb. Use of Vt made 
no detectable difference in overall patterns. Some minor increases 
in value in the vicinity of sharp curvature at low levels was noted but 
were of the order of only 10 to 20 percent. Advection of vorticity with 
w’xwas not computed as its significance has been seen to be small by 
its small effect when added to Ve. However, advection of temperature 
using Vx was computed using a program which evaluates the expressions 



J (E,T) 


(21) 


\V *7 T . 


(22) 



One hour change of temperature was computed. Figure 14 shows 1000= 
mb change of temperature with Ve (solid lines) . Units are degrees 
Celsius lO”^. Detailed analysis of the grid map shows that advection 
with Vx over the entire map was less than ^ 0.01 degrees over the United 
States except in the area along the Gulf Coast where values are depict- 
ed with a dashed contour enclosing an area of - 0.2 degrees. Here the 
advection with Vx actually exceeded EADV. This is the area noted earl- 
ier of the deepening low pressure system and front. Relatively large 
values were recorded at the 850-mb and 700-mb levels in this same area. 
Again the rapidity of movement of the front and associated values of ad- 
vection and short length of the data period prevented detailed investiga- 
tion but suggests the possibility of a meteorologically significant rela- 
tionship. 
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5. Conclusions 



1. Where there is curvature in the pressure patterns involved, the non 
divergent velocity is a better approximation to the wind than the geo- 
strophic value. 

2. No significant change is made by combining the non-divergent and 
divergent winds for the uses investigated. 

3. A relationship possibly exists between the divergent wind and front 
al areas that could be utilized in locating such phenomena. This con- 
clusion is based on very limited data but offers an interesting avenue 
for further investigation. 



15 



FIGURE 1 




16 



FIGURE 2 
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APPENDIX I 



Pillow is defined as 

X'-h = 3^68 







(I-l) 



The root-mean- square differnece (RMSD) is determined as follows: 



x*n=3e(i,g ^ 

. (l-^) 

In equations (I-l) and (1-2), and represent the point values of 
the fields being compared. In this example the A value is a geostrophic 
value of speed or direction and B is the corresponding value of either non- 
divergent or total wind. 
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